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A large fraction of the information collected by cosmological surveys is simply discarded to avoid 
lengthscales which are difScult to model theoretically. We introduce a new technique which enables 
the extraction of useful inforination from the bispectrum of galaxies well beyond the conventional 
limits of perturbation theory. Our results strongly suggest that this method increases the range 
of scales where the relation between the bispectrum and power spectrum in tree-level perturbation 
theory may be applied, from fcmax ~ OAhMpc~^ to ^ 0.7/iMpc~^. This leads to correspondingly 
large improvements in the determination of galaxy bias. Since the clipped matter power spectrum 
closely follows the linear power spectrum, there is the potential to use this technique to probe the 
growth rate of linear perturbations and confront theories of modified gravity with observation. 



Introduction. — At the largest observable scales, den- 
sity fluctuations in the Universe have sufRciently low 
amplitude to allow a simple approximation to model 
their evolution, in the form of linear perturbation the- 
ory. But at progressively smaller scales this formalism 
breaks down, and more complex prescriptions are re- 
quired to predict the clustering patterns. While in re- 
cent years there has been steady progress improving our 
understanding of the nonlinear matter power spectrum 
(e.g. [1-^]), such developments are unlikely to keep pace 
with the precision demanded by the ambitious galaxy 
surveys of the future. If we are to take full advantage of 
their capabilities, these surveys will require knowledge of 
the matter power spectrum to sub-percent precision at 
lengthscales well into the nonlinear regime. 

The galaxy bispectrum is able to offer insight into key 
topics in cosmology. Along with a theoretical model of 
its evolution, a measurement of the bispectrum enables 
properties of the galaxy bias to be determined [4]. It 
also potentially holds a signature of inflation, in the form 
of primordial non-Gaussianity. However, like the power 
spectrum, the bispectrum is known to be strongly in- 
fluenced by nonlinearities. To obtain robust parameter 
constraints, it is usual to discard information above a lim- 
iting /cmax, an obviously undesirable state of affairs. At 
low redshifts we are restricted to working on extremely 
large scales, typically fcmax ^ 0.1 /iMpc~^ [5[, where cos- 
mic variance severely hampers our precision. This value 
of fcniax niay be improved with more detailed modelling, 
going beyond tree-level perturbation theory to one-lo_o 
or Renormalised Perturbation Theory [3|, |5| ; see also 
These methods are able to extend the /c-range of ap- 
plicability by a moderate amount over tree-level, at the 
cost of considerable complexity. An alternative to mod- 
elling the full matter power spectrum is to manipulate the 
real space density field such that the effects of nonlinear 
growth are suppressed, an approach that has previously 



been demonstrated to successfully enhance the baryon 
acoustic oscillations (BAO) [7|. There is, therefore, mo- 
tivation to study such transformations as a general tool 
for cosmological analysis. 

The simplest class of reconstruction methods are lo- 
cal, monotonic mappings from the initial to final density 
fields. The Gaussianisation process proposed by Wein- 
berg [8| enforces Gaussianity in the field's one-point dis- 
tribution, and this has been shown to recover the shape of 
the linear power spectrum [9|]. Some other filters, includ- 
ing log(l -f 6) [10|, perform comparably well. What such 
transformations have in common is the penalisation of 
the highest density regions of the field, a feature we shall 
exploit further. Our goal is to extend the range of the 
relatively simple tree-level perturbation theory by ma- 
nipulating the field in real space. This work begins with 
the ansatz that the nonlinear behaviour of the density 
field is localised not only at high wavenumbers in Fourier 
space, but also within high density regions in real space. 

Peak Clipping. — Extreme peaks of the density field 
correspond to the regions furthest from the quasi-linear 
evolution we wish to study, yet it is these peaks that dom- 
inate the contribution to higher-order statistics. There- 
fore, a promising approach for extracting the primordial 
or tree-level gravitational bispectrum is to 'clip' the den- 
sity field, reducing the most extreme peaks within a sam- 
ple. We maintain a continuous density field by reducing 
the value of extreme density peaks down to a threshold 
value, such that we generate the truncated field 6t{x) 



Stix)^Six), (<5(x) < J--'^) 
^t(a;)=<5'"^'^, {6{x) > J'"^'^) 



(1) 



where the density contrast is defined by S{x) = 
p{x)/p{x) — 1. We then re-evaluate St{x) using the trun- 
cated mean density to ensure {St{x)) — 0. This pre- 
scription allows us to generate results which are rather 



insensitive to the choice of threshold S"^'^^, provided that 
it is low enough to affect at least ~ 0.1% of the volume. 

Galaxy redshift surveyors are well accustomed to deal- 
ing with effects that arise from the process of astronomi- 
cal observation. For measurements of statistics in Fourier 
space, these are conventionally eliminated by deconvolu- 
tion of a window function — yet here the effective mask 
cannot be treated in this manner, as it is strongly corre- 
lated with the underlying density field. Crucially, how- 
ever, the mask covers such a small fraction of the survey 
volume that we do not need to correct for it in this work. 

Bispectrum. — We begin with the assumption that on 
large scales the number density of galaxies is well mod- 
elled in terms of the local linear and quadratic bias pa- 
rameters 61 and 62, defined such that Sg{x) = biS{x) + 
^S^{x), where S and 6g denote the dark matter and 
galaxy fields respectively. 

In a given sample of halos or galaxies, we can attempt 
to recover the bias properties by comparing the observed 
bispectrum Bg with the tree-level theoretical expectation 
Bg. These two functions are (see e.g. |lll4l4| ) 



for the window function imposed by the NGP method 
(see e.g. [l9| ) 



sik)=5o{k) n 
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2 sin 






(5) 



where fcnyq is the Nyquist frequency, and p = 1 for NGP. 
The fundamental frequency of the box {kf ^ 
0.013 ft.Mpc~ ) dictates the bin width for estimating our 
statistics. We perform a brute-force evaluation of all tri- 
angle configurations up to kmax — T^^kf ~ 0.9/iMpc^ , 
storing the cumulative values as a function of the vector 
amplitudes rounded to the nearest multiple of kf . 



C{ka, kt, k,) = J2Y. S{ka)5{h)S{k, 



ka fcb 



(6) 



where kc = — ika + kbj . This raw cumulative count is 
rapidly distilled to any desired configuration 



{5g{kl)Sg{k2)Sg{k3)) = (271)^5^ (fcl , ^2 , ^3 )'5(^1 + ^2 + ^3 ) 

(2) 
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2j{h,k,) + ^ Pg{h)Pgik,), 

(3) 

where P{k) = {\6{k)\'^), and the galaxy power spectrum 
is Pg{k) = b\P{k). J{ki,k2) is given by 



J{ki,k2) 



1-D{nrn) + ^- [-jT + jr] +D{nrn) COS^ ' 



1 _ 3 0-1/143 

2 14 "m 



where 9 denotes the angle between fci and ^2. This has 
a weak dependence on cosmology. For a flat ACDM Uni 
verse, D may be modelled by D{V,m) 

^m. 

Numerical Methods. — We use data from the 
Millennium-I simulation, utilising both the dark matter 
density field and an associated catalogue of ~ 7 million 
galaxies [17J, |l8|. The redshift zero snapshot presents a 
strongly evolved density field (with as = 0.9), providing 
a suitably challenging test for this method; galaxies are 
selected by the mass cuts logio(M»/M0ft.) > 9 or > 10. 
These cuts yield number densities of 5.8 x 10~^ h^Mpc^'^ 
and 1.3 x 10~^ /i^^Mpc"^ respectively. 

A galaxy number density field is constructed in the 
same manner as the dark matter density field, at the 
same resolution (256'^; box size — 500/i~^Mpc) and by 
the Nearest Grid Point (NGP) algorithm. This field is fil- 
tered in real space in accordance with ([T]) , with S'^'^^ cho- 
sen such that 0.1% of the volume is clipped. We correct 



where we define ki > k2 > k^, and nt is the total number 
of triangles counted. This is evaluated alongside P{ki), 
P{k2), Piks) and the mean configuration J{ki,k2)- We 
select a tolerance value such that the vector magnitudes 
were an appropriate proportion, matching the desired 
values of ki/k2 and ki/k^ to within 10%. We assume 
diagonal covariance of the bispectrum estimates, deter- 
mined by jackknife resampling from eight subsamples. 
This is a potentially important issue which we will ex- 
plore in future work, as it could result in underestimation 
of the error bars, but the evidence from the bias recovery 
in Figure [2] is that it isn't a significant effect in this study 

In order to determine whether this technique has suc- 
cessfully recovered the properties of the linearised field, 
we need to estimate the true values of 6i and 62. For the 
dark matter, this is simply given by 61 = 1, 62 == 0. For 
the galaxy field, we fix 5i from the ratio of large-scale 
power spectra, and determine 62 by least-squares fitting 
to fields smoothed with a cutoff fcmax = 0.2 ft,Mpc~^. 

For a galaxy field 61 and 62 are not invariant under the 
clipping process. Galaxy bias is sensitive to the selection 
criteria of the sample, and so we find that for truncation 
fractions around ^1% the underlying bias values change 
by up to ~ 10%. Note that the change of bias is not 
necessarily important: a useful by-product of this proce- 
dure is to determine the matter power spectrum, which 
can be obtained from the clipped galaxy power spectrum 
and the clipped bias. An advantage is that the resulting 
clipped matter power spectrum is very close to the linear 
power spectrum. 
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FIG. 1. Left: The left and middle panels show the dark matter bispectrum ([2)) for colinear triangles with wave vectors in the 
ratio 3:2:1 and 2:1:1 (left and middle). The solid lines are estimates of the bispectrum using the matter power spectrum, as 
given by (|3]). The upper set of points and line relates to the untouched density field, while the lower points and lines corresponds 
to a field which has first been clipped in accordance with ^. The truncation limits are chosen such that 0.1% (middle) and 1% 
(lower) of the volume is clipped. Right: The same as the left-hand panel, but here we study the logiQ(Mt/M0/i) > 9 galaxy 
sample. 



Results. — Figure [T] displays the impact of clipping on 
the bispectrum of the dark matter density field (left and 
middle panels) and the galaxy number density field (right 
panel). In each panel, the upper set of data points rep- 
resent the undipped bispectrum, while the middle and 
lower sets of data points show the bispectrum after 0.1% 
and 1% of the field has been clipped. The solid lines 
denote the tree-level prediction for the bispectrum. The 
effect of the clipping is to reduce the extreme nonlineari- 
ties, bringing the tree-level prediction into much closer 
agreement with the data and permitting much higher 
wavenumbers to be incorporated. The left and right pan- 
els have wave vectors in the ratio 3:2:1, and the middle 
panel 2:1:1. Beyond k > 0.7 hMpc"^ the correction given 
by ([S]) becomes large and the results in this regime may 
not be reliable. 

Clipping 0.1% of the volume corresponds to a threshold 
oiSmax — 60 and 70 for the dark matter and galaxy fields, 
corresponding to mass fractions of ^ 14% (dark matter) 
and ~ 15% (galaxies). As the fraction of the volume 
subject to clipping increases, so does the shot noise, as 
the number density of galaxies and the amplitude of the 
filtered power spectrum are both suppressed. Note that 
coarser grid spacings would require larger fractions of 
the volume to be clipped, in order to ensure a sufficient 
number of extreme peaks are affected. 

Any given configuration i?(fci, A;2, ^3) exhibits degen- 
eracy between the parameters 61 and &2 ■ Figure [2] illus- 
trates how a combination of two configurations (3:2:1 and 
3:2:2) helps to lift this degeneracy. We compute the like- 
lihood of 61 and 62 from equation ^, assuming that the 
errors in the power spectrum are negligible in compari- 
son with the bispectrum. The red and blue contours are 



derived from undipped and clipped (0.1%) density fields 
respectively, and represent the one- and two- a confidence 
contours. The clipping allows us to reach a higher fcmax, 
substantially reducing the errors on the bias parameters. 
The left panels correspond to the dark matter field, while 
the middle and right panels are for galaxy fields with 
mass cuts \ogiQ{M^ / MqH) > 9 and > 10 respectively. 

For the original undipped field we take fcmax — 
0.17/iMpc~^. If we attempt to enlarge this value, the 
confidence contours rapidly diverge from the true value. 
When working with the filtered field, we are able to utilise 
the bispectrum up to 0.5 hMpc^^ for the scalene flattened 
(3:2:1) shape. The isosceles (3:2:2) is slightly less respon- 
sive to the clipping process, and for that configuration we 
apply a cut at Amax = 0.35 hMpc'^. We note that higher 
truncation fractions appear to allow even higher values 
of fciiiax for the co-linear configurations. Triangle con- 
figurations closer to the equilateral shape suffer from a 
decrement at larger truncation fractions. We speculate 
this may be due to the co-linear Fourier triangles pick- 
ing out pancake-like structures, which are lower density 
and thus less susceptible to clipping than the filaments, 
to which equilateral configurations are more sensitive. It 
is expected that this effect could be calibrated using N- 
body simulations. 

We have also repeated these tests with a higher red- 
shift {z = 0.687) dark matter field, effectively probing 
the dependence on the power spectrum. We find similar 
results for these cases, that clipping just 0.1% of the field 
significantly extends the reach of (|3]). We investigated 
lognormal mapping and Gaussianisation, which can be 
effective at power spectrum level [8|-[l0j; however at bis- 
pectrum level we did not find good agreement with ([3]). 




FIG. 2. Constraints on the linear and quadratic bias parameters, as derived from combining two configurations of tfie bis- 
pectrum, before (red) and after (blue) clipping 0.1% of the volume. Left: This is for the z = dark matter density field, 
so the correct values are 6i = 1 and 62 = 0, as denoted by a black cross. For the flat configuration, the clipping allows us 
to approximately triple the Fourier cut fcmax from 0.17 /iMpc"'^ to 0.5 /iMpc~^. The dotted contours outline the individual 
contributions for the case of the filtered field. Middle and Right: Same as the left panel, except here we are using a mock 
galaxy catalogue with mass cuts logjQ(M»/M0/i) > 9 and > 10 respectively. 



Conclusions. — Determining the form of the dark mat- 
ter power spectrum remains a principal goal of modern 
cosmology. We have demonstrated how the simple pre- 
scription of clipping high density peaks allows a much 
larger volume of fc-space to be incorporated into such 
predictions. This leads to a substantially more accu- 
rate determination of the bias, and consequently tighter 
constraints on cosmological parameters. For the highly 
evolved field of the Millennium simulation (erg — 0.9 at 
redshift zero), we find that the maximum wavenumber 
fcmax may typically be extended by this technique from 
-- 0.1 /iMpc"^ to beyond - 0.5 /iMpc"\ The fuh impact 
of this extension upon signal to noise becomes apparent 
upon noting that the number of triangles contributing 
to the bispectrum computations scales in proportion to 

'*'inax' 

In this letter we have chosen to focus on the galaxy bias 
parameters, but it should be clear that this technique is 
extensible to quite generic measurements in cosmology. 
At the heightened precision with which the matter bis- 
pectrum can now be studied, it is natural to consider 
whether improved constraints on f^^ sue attainable. Al- 
though any signature of primordial non-Gaussianity is 
expected to be most prominent at the largest scales, the 
abundance of modes at high k may compensate for the 
smaller signal. The clipping procedure may also prove 
valuable in improving cosmological constraints on the 
neutrino mass [20i - i2S] . and in modified gravity models 
where the linear growth rate differs from General Rela- 
tivity. 

There are a number of ways in which the technique 



could be developed further. A closer-to-optimal approach 
is to evaluate the relationship between the preferred trun- 
cation threshold and particular wavenumber values. In- 
deed the ultimate limit on fcmax is at present uncertain, 
since larger truncation fractions extend the range of va- 
lidity beyond 0.7 hMpc~ . The compensation of the clip- 
ping window, and moving to redshift space, where there 
wiU undoubtedly be challenges due to 'Fingers of God', 
will be topics of forthcoming work. Exploitation of the 
bispectrum to the level of precision that will soon be ob- 
servationally accessible will likely require the incorpora- 
tion of scale-dependence into the bias model. Considering 
these future directions, it seems reasonable to conclude 
that density filtering schemes may establish themselves 
as an important tool for cosmological analysis. 
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